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ABSTRACT 

A relaxation method is demonstrated which reliably solves the nonlinear two-point-boundary-value 
problem which arises when optimal control theory is applied to determination of large angle 
maneuvers of flexible spacecraft. The basic ideas are summarized and several idealized maneuvers 
are determined. The emphasis is upon demonstrating the basic ideas and practical aspects of the 
methodology. References are cited, particularly Turner’s dissertation which presents detailed 
formulations and more general applications. 
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Discussion of Figures 


With reference to Figure 1 , we employ the method of assumed modes to obtain a set of ordinary 
differential equations which govern deflections and rotations. The form of the equations of motion 
are given in Figure 2. Note the high dimensionality and the variability of the coefficient matrix. 
Note that solution for the acceleration coordinates is required in order to integrate motion as a 
function of time, and in order to apply optimal control theory. 

Figure 3 displays a partitioned algorithm which efficiently determines the inverse of the high- 
dimensioned, configuration-variable coefficient matrix. Consistent with this partitioning algorithm, 
we consider in Figure 4 an algorithm for obtaining partial derivatives of the inverted coefficient 
matrix with respect to deflection coordinates (required in the optimal control algorithm). 

Figure 5 summarizes the state and co-state differential equations which follow from Pontryagin’s 
principle as the necessary conditions satisfied by optimal (minimum quadratic cost) maneuvers. 
Observe that the initial and final states are generally known, but the initial and final co-states are 
usually unknown. Thus, as usual, a nonlinear two point boundary value problem (TPBVP) has 
resulted. Notice the quadratic angular velocity nonlinearity due to “rotational stiffness.” 

In Figure 6, we summarize an imbedding/relaxation approach which has proven a reliable approach 
for solving TPBVP’s of the above structure. In essence, a one parameter (a) family of problems is 
constructed that one special member (cv = 0) has an analytical solution, while another member 
(a = 1) is the true problem of interest. By relaxing oc through a sequence of increasing values 
0 < 1, we can extrapolate arbitrarily good initial or final co-state estimates (by adjusting the 

a-increment) from previous converged solutions, thereby allowing efficient differential corrections 
to isolate accurate co-states corresponding to each a. Typically, only 4 or 5 ai values are actually 
required to reach the desired oc = 1 solution. This method and related methods are developed and 
applied to several examples in Reference 3. 

Considering now a specific configuration, we refer to Figure 7. The four identical cantilevered 
appendages are mounted in the same plane to the rigid central hub. We neglect the hub radius 
in any equation in which it appears divided by the appendage length. Referring to Figure 8, we 
restrict attention to pure spin rotations and antisymmetric deflections, consistent with spin-up, 
spin-down, and rest-to-rest maneuvers with the configuration initially and finally undeformed. We 
consider only the case of torques applied to the hub. 

Table 1 describes seven maneuver calculations, corresponding to three sets of maneuver boundary 
conditions and four different dynamical models. These cases are selected to demonstrate the 
effects of rotational stiffening and to show that the relaxation method can handle both high 
dimensionality and nonlinearities. 

Figures 9a - c display the angle of rotation, angular rate and torque for the case 1 maneuver (rigid 
appendages). For comparison, Figures 10a - c display the same variables for cases 2L and 2N of 
flexible appendages, assuming a 1 mode expansion. It is of interest to note that the flexibility 
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effects are large indeed. The flexible case torque oscillates anti-symmetrically about the rigid case 
torque, the desired final angle and angular rate are achieved and the modal amplitude (and its de- 
rivative) are simultaneously driven to zero. It is interesting that the linear and nonlinear solutions 
were identical, to graphical accuracy, due to the small deflections and velocities of this particular 
maneuver. 

Figure 1 la - d and 12a - d display angle of rotation, torque history, and amplitudes of the first two 
modes for cases 3L and 3N, respectively. The maneuver is an extremely rapid spinup from rest to 

0.5 rad/sec in 60 sec. The linear (3L) and nonlinear (3N) solutions differ significantly, but the 
linear solution retains the general shape and amplitudes differ by less than 10% throughout most of 
the motion. 

Figure 13a - g display the angle of rotation, angular rate, torque, and the first four modal amplitudes 
for case 4L (a rest-to-rest maneuver through a 360"^ rotation). These results simply show that, 
indeed, the large rigid rotations and vibration suppression of several degrees of freedom are deter- 
mined. 

We offer the following significant conclusions: 

• An Optimal Control Formulation is Presented for General 3 Dimensional Maneuvers of a 
Class of Flexible Satellites 

• A Partitioning Method is Introduced to Invert the Ro tati on al-Vi bra ti on al Equations of 
Motion for Acceleration Coordinates and to Obtain the Adjoint Equations 

• An Imbedding/Relaxation Process if Demonstrated for Solution of the Two-Point-Boundary- 
Value Problem. 

• Numerical Studies Indicate that Practical Algorithms Result from these Developments 
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THE METHOD OF ASSUMED MODES 


The deflection of the jth flexible member is modeled as 


Uj(x,y,z,t) = 

i = l 
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The sets of spatial "assumed modes" 

|u^^(x,y,z) ... . (x,y,z) j- . . .| U^^(x,y,z) ... U^.(x,y,z)| 

{ V^^(x,y,z) ... V^^.(x,y,z) V^^(x,y,z) ... V^^.(x,y,z)| 

^ W-] 1 (x»y»z) ... W.|^(x,y,z) |...<| W^.|(x,y,z) ... W^^{x,y,z)| 

are prescribed . As minimum requirements, they must 

• be linearly independent 

• satisfy u, v, w's geometric boundary conditions 

The amplitude functions constitute the aonf^lgu/iatyLon vtatoH. 
in(^) 



The amplitude's play the role of discrete generalized coor 
di nates. 

Figure 1 
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ROTATIOfWyDEFLECTIONAL EONS OF MOTION 

Dynamics of flexible satellites 


{0} = [F(0)]{u)} 

+ {f (0»(»),r|>ri»t)} + {u} 


( 1 ) 

( 2 ) 


[M]{r)} * + {g(6,w,ri,ri j t) } 


(3) 


Combine (2) & (3) 


( J(n) 
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( 4 ) 


Note 


[H] & [M] are constant 

[J(n)] = [J^] + [Jjj(n)] , l|jp|| < < IIJqII (typically) 

Inertia of Inertia varia- 

unde formed tlons due to 

vehicle deformations 


A problem: 

We need eqns of motion in the state space form x = F(x,u,t), 
but 


(i) The coefficient matrix of (4) is variable 
(ii) Its dimensions may be several hundred 


Figure 2 
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PARTITI0I€D/PERHJI«AT1CIN INVERSION OF 1TC 
COEFFICIENT mTRlX 


Name the submatrices: 


^11 

^ 2 } 

^21 

^22 




The C.. can be expressed directly as a function of J, H, H as: 

* w 



Form 1 

Form 2 


C 2 ^ 

T -1 

(0 - h‘ M ' H) 

C22 

1 T -■* 

(M - H J ' h') 

M’^ - M"^ H 02 ! 

^21 

-C 22 H 0-' 
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For direct numerical calculations. Form 2 is preferred since 
(i) (J - M“^ H) is a 3 X 3 matrix 

(ii) M is generally diagonally dominant (an identity matrix if one first 
solves an eigenvalue problem - Note M is positive-definite symmetric) 


Figure 3 



CERTAIN REQUIRED PARTIAL DERIVATIVES S HOW TO DETERMINE THEM 

Rotational/Vibrational Equations of Motion 
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Figure 4 
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FO(flUL'\TION OF THE OPTICAL COTOL PROBLB1 


STATE VARIABLES 

*2 “ —2 “ ^4 “ 

STATE DIFFERENTIAL EQUATIONS 
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Find ^(t) generating a trajectory initiating at terminating at 

which minimizes the function 


J - I (»■' H j 4 “11 il'"' 

o i=2 


HAMILTONIAN 


H = i=2 j, 4 


4 

I 
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PONTRYAGIN'S NECESSARY CONDITIONS 


Co-State Equations 

aH _ 




X - - 


Optimal Control 

Minimize H at each instant with respect to admissible u(t), this 
yields — U > • • > y X ^ ^ * • y X^ ^ t) 


Figure 5 
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I^EEH)INI5/REUWATI0N METHOD FOR SOLVING 
TWO-POINT BOUNDARY VAIiJE PRffiliM 


Define merged vector 



The coupled state and costate differential equations are then 
Z_= [A]Z^ + a{all nonlinear terms} 


Typically, we know x(tQ) and x(t^), but not i(tQ), x(t^). 

For a= 0, we can solve for X.(tQ) exactly. 

By taking sufficiently small a-increments, we can use converged lit^) from 

neighboring optimal solutions to initiate successive approximations with 
oAbitAjvUly good -itoAdxng e^timate^ for the unknown A(tp). 

Typically, only 5 to 10 intermediate a-values are required a practical 
algorithm results. 


Figure 6 
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Figure 9 Case 1 Rigid Body Figure 10 Cases 2L, 2N Flexible 
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Figure 13 Case 4L Rest—To—Rest Maneuver 
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TABLE 1 DESCRIPTION OF TEST CASE MANEUVERS 


Case # 

Qualitative Description 

of Modes (N) 
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1 

Rigid Appendages 
Rest-to-Rest Maneuver 
t^ = IA.221 sec. 

0 
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2L 

Linear Kinematics 
Rest-to-Rest Maneuver 
t^ = 2ir/u>j^ = 14.221 sec 
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2N 

Nonlinear Kinematics 
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Linear Kinematics 
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Nonlinear Kinematics 
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